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Summary 


Swimming  efficiencies  of  fish  and  cetaceans  have  been 
related  to  a  certain  synchrony  between  stroke  cycle 
frequency,  peak-to-peak  tail/fluke  amplitude  and  mean 
swimming  speed.  These  kinematic  parameters  form  a  non- 
dimensional  wake  parameter,  referred  to  as  a  Strouhal 
number,  which  for  the  range  between  0.20  and  0.40  has 
been  associated  with  enhanced  swimming  efficiency  for 
fish  and  cetaceans.  Yet  to  date  there  has  been  no  direct 
experimental  substantiation  of  what  Strouhal  numbers  are 
preferred  by  swimming  cetaceans.  To  address  this  lack  of 
data,  a  total  of  248  Strouhal  numbers  were  calculated  for 
the  captive  odontocete  cetaceans  Tursiops  truncatus, 
Pseudorca  crassidens,  Orcinus  area,  Globicephala  melaena, 
Lagenorhynchus  obliquidens  and  Stenella  frontalis. 
Although  the  average  Strouhal  number  calculated  for 
each  species  is  within  the  accepted  range,  considerable 
scatter  is  found  in  the  data  both  within  species  and  among 
individuals.  A  greater  proportion  of  Strouhal  values  occur 
between  0.20  and  0.30  (74%)  than  the  0.25-0.35  (55%) 
range  predicted  for  maximum  swimming  efficiency. 
Within  0.05  Strouhal  increments,  the  greatest  number  of 
Strouhal  values  was  found  between  0.225  and  0.275  (44%). 


Where  propulsive  efficiency  data  were  available  (Tursiops 
truncatus,  Pseudorca  crassidens,  Orcinus  area),  peak 
swimming  efficiency  corresponded  to  this  same  Strouhal 
range.  The  odontocete  cetacean  data  show  that,  besides 
being  generally  limited  to  a  range  of  Strouhal  numbers 
between  0.20  and  0.40,  the  kinematic  parameters 
comprising  the  Strouhal  number  provide  additional 
constraints.  Fluke-beat  frequency  normalized  by  the  ratio 
of  swimming  speed  to  body  length  was  generally  restricted 
from  1  to  2,  whereas  peak-to-peak  fiuke  amplitude 
normalized  by  body  length  occurred  predominantly 
between  0.15  and  0.25.  The  results  indicate  that  the 
kinematics  of  the  propulsive  fiukes  of  odontocete  cetaceans 
are  not  solely  dependent  on  Strouhal  number,  and  the 
Strouhal  number  range  for  odontocete  cetaceans  occurs  at 
slightly  (~20%)  lower  values  than  previously  predicted  for 
maximum  swimming  efficiency. 

Key  words:  Strouhal  number,  swimming,  dolphin,  cetacean, 
odontocete,  Tursiops,  Pseudorca,  Orcinus,  Lagenorhynchus, 
Globicephala,  Stenella,  Delphinapterus. 


Introduction 

An  important  aspect  of  swimming  is  the  ability  to  move 
efficiently.  Paradoxically,  early  attempts  at  building  fish- 
inspired  mechanisms  achieved  disappointingly  low 
efficiencies  (Triantafyllou  and  Triantafyllou,  1995).  It  was 
only  through  a  deeper  understanding  of  the  vorticity  produced 
along  the  swimming  animal  and  in  its  wake  that  significant 
progress  was  achieved.  Beginning  almost  40  years  ago,  Rosen 
(1959,  1961,  1963)  discerned,  through  a  series  of  innovative 
flow  visualization  experiments,  a  system  of  vortices  appearing 
along  the  sides  of  swimming  fish  and  dolphins.  Rosen  (1959, 
1961,  1963)  hypothesized  that  some  of  the  rotational  energy 
surrounding  the  undulating  motion  of  a  fish  or  dolphin  could 
be  regained  for  propulsion  through  proper  synchronization  of 
the  animal’s  body  to  the  vortex  flow.  Rosen  (1959,  1961) 
further  deduced  an  equation  for  fish  and  dolphin  motion.  This 
equation  predicted  swimming  speed  to  be  proportional  to  the 
product  of  the  tail  beat  amplitude  and  frequency.  Rosen  (1959, 


1961)  referred  to  this  proportionality  as  the  ‘fish’  constant  and 
hypothesized  that  it  was  nearly  the  same  for  fish  and  dolphins. 

A  similar  conclusion,  but  through  more  rigorous  theoretical 
analysis  and  detailed  experimental  studies,  has  been  drawn  by 
Triantafyllou  et  al.  (1991,  1993).  Performing  stability  analysis 
of  the  mean  velocity  profiles  of  a  pitching  airfoil,  Triantafyllou 
et  al.  (1991,  1993)  have  shown  that  maximum  spatial 
amplification  and  optimum  creation  of  thrust-producing  jet 
vortices  lies  in  a  narrow  range  of  nondimensional  frequencies 
referred  to  as  the  Strouhal  number  (Sf).  The  predicted  St  range 
for  maximum  spatial  amplification  occurs  between  0.25  and 
0.35,  peaking  at  0.30  (Triantafyllou  et  al.,  1991,  1993; 
Streitlien  and  Triantafyllou,  1998).  Triantafyllou  and 
Triantafyllou  (1995)  have  argued  that  for  5f=0.25-0.35, 
swimming  efficiency  for  fish  and  cetaceans  would  also  peak. 
Experiments  with  isolated  oscillating  foils  have  found  highest 
propulsive  efficiencies  for  St  between  0.20  and  0.40 
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(Triantafyllou  et  al.,  1991,  1993,  2000;  Anderson  et  al.,  1998; 
Read  et  al.,  2003). 

The  St,  which  is  related  to  how  fast  the  vortices  are  being 
generated  and  the  space  between  them,  is  defined  as: 

St  =  Af/U,  (1) 

where  A  is  the  width  of  the  wake,  taken  to  be  equal  to  the  peak- 
to-peak  maximum  excursion  of  the  trailing  edge  of  the  foil  or 
fluke,  /  is  the  frequency  of  oscillation,  and  U  is  the  mean 
forward  velocity.  St  is  in  fact  the  inverse  of  the  ‘fish’  constant. 
Similar  combinations  of  these  kinematic  swimming  parameters 
have  previously  been  made  to  characterize  the  swimming 
motion  of  fish  (e.g.  Rosen,  1959,  1963;  Pyatetskiy,  1970; 
Webb,  1975),  cetaceans  (Semonov  et  al.,  1974;  Kayan  and 
Pyatetskiy,  1977)  and  athletes  (Pershin,  1988),  but  not  within 
such  a  predictive  theoretical  framework. 

Although  many  species  of  cetaceans  are  believed  to  be 
particularly  swift,  efficient  swimmers,  corresponding  St  data 
have  been  extremely  limited.  The  cetacean  Strouhal  number 
data  of  Triantafyllou  et  al.  (1993)  consist  of  only  two  values. 
Moreover,  they  were  indirectly  derived  from  analysis  of  traces 
(Lang  and  Daybell,  1963)  obtained  from  motion  picture  frames 
of  a  2.03-m  Pacific  white-sided  dolphin  {Lagenorhynchus 
obliquidens)  swimming  non-uniformly  within  a  tank.  One  St 
value,  0.32,  corresponded  to  the  dolphin  swimming  while 
wearing  a  1.91  cm-diameter  drag  collar.  The  remaining  St  of 
0.30  corresponded  to  swimming  without  the  drag  collar. 
References  to  this  less  than  optimal  data  set  (Triantafyllou  et 
al.,  1991,  1993,  2000;  Triantafyllou  and  Triantafyllou,  1995; 
Taylor  et  al.,  2003)  appear  repeatedly  throughout  the  literature. 

Many  different  species  of  captive  odontocete  cetaceans  have 
been  trained  to  swim  steadily  behind  viewing  panels  and 
provide  a  unique  opportunity  for  a  much  larger,  more  accurate 
St  database.  Here,  Strouhal  numbers  are  calculated  from 
recordings  of  six  species  of  trained  odontocete  cetaceans.  The 
species  include  members  of  the  family  Delphinidae:  bottlenose 
dolphin  (Tursiops  truncatus),  false  killer  whale  (Pseudorca 
crassidens),  spotted  dolphin  (Stenella  frontalis),  striped 
dolphin  {Lagenorhynchus  obliquidens),  killer  whale  {Orcinus 
orca)  and  pilot  whale  {Globicephala  melaena).  For 
comparison,  Strouhal  numbers  for  the  slower,  less  efficient 
swimming  beluga  whale  {Delphinapterus  leucas',  family 
Monodontidae)  are  also  included  but  treated  separately. 

The  objective  of  the  present  study  was  to  investigate  the 
range  of  Strouhal  numbers  preferred  by  swimming  cetaceans, 
how  this  range  varied  between  species,  within  species  and  for 
individuals  and,  most  critical,  what  range  of  Strouhal  numbers 
corresponded  to  maximum  propulsive  efficiency.  Data  directly 
relating  swimming  efficiency  to  St,  for  any  swimming  animal, 
have  previously  not  existed  (Bandyopadhyay  et  al.,  2000). 

Materials  and  methods 

Experimental  animals 

The  swimming  motions  of  seven  species  of  trained 
odontocete  whales  were  recorded  at  Sea  World  in  Orlando,  FL, 


USA,  San  Antonio,  TX,  USA  and  San  Diego,  CA,  USA  and 
the  National  Aquarium  in  Baltimore,  MD,  USA.  The  species 
examined  were  from  the  family  Delphinidae,  including  one 
spotted  dolphin  {Stenella  frontalis  Cuvier),  two  Pacific  white¬ 
sided  dolphins  {Lagenorhynchus  obliquidens  Gill),  11 
bottlenose  dolphins  {Tursiops  truncatus  Montagu),  four  false 
killer  whales  {Pseudorca  crassidens  Owen),  one  pilot  whale 
{Globicephala  melaena  Lesson)  and  six  killer  whales  {Orcinus 
orca  Linnaeus),  and  the  family  Monodontidae,  including  three 
beluga  whales  {Delphinapterus  leucas  Pallus).  Some  of  the 
data  were  previously  reported  in  Fish  (1993,  1998). 
Morphological  measurements  for  each  animal  can  be  found  in 
Table  1.  Body  length,  L,  is  defined  as  the  linear  distance  from 
the  rostral  tip  to  the  fluke  notch.  The  animals  subsisted  on  a 
diet  of  herring,  smelt,  mackerel  and  squid  supplemented  with 
vitamins,  dispersed  at  irregular  intervals  throughout  the  day. 

Experiments  were  performed  in  large  elliptical  pools  with 
maximum  lengths  of  27.4^8.8  m.  The  curved  portions  of  each 
pool  were  constructed  of  1. 7-2.1  m-wide  Plexiglas  panels 
separated  by  0.2  m-wide  posts,  allowing  for  an  unobstructed 
view  of  the  animals  as  they  swam.  A  water  depth  of  1.4— 2.1  m 
was  visible  through  the  panels.  The  depth  of  the  pools  was 


Table  1.  Morphometries  of  odontocete  cetaceans 


Species 

Sex 

Body 

mass 

(kg) 

Body 

length 

(m) 

Age 

(years) 

Stenella  frontalis 

Male 

81.6 

1.83 

5 

Lagenorhynchus 

Male 

130.6 

2.18 

16 

obliquidens 

Lagenorhynchus 

Female 

139.7 

2.25 

15 

obliquidens 

Tursiops  truncatus 

Female 

146.9 

2.49 

4 

Tursiops  truncatus 

Male 

149.2 

2.06 

3.5 

Tursiops  truncatus 

Female 

187.8 

2.57 

17 

Tursiops  truncatus 

Female 

192.8 

2.51 

20 

Tursiops  truncatus 

Male 

211.4 

2.54 

18 

Tursiops  truncatus 

Female 

212.3 

2.59 

8.5 

Tursiops  truncatus 

Female 

219.1 

2.69 

16 

Tursiops  truncatus 

Male 

228.2 

2.70 

20 

Tursiops  truncatus 

Female 

229.5 

2.60 

20 

Tursiops  truncatus 

Female 

256.3 

2.95 

13 

Tursiops  truncatus 

Female 

263.1 

2.70 

22 

Pseudorca  crassidens 

Female 

358.3 

3.57 

9 

Pseudorca  crassidens 

Female 

391.3 

3.71 

10 

Pseudorca  crassidens 

Female 

508.1 

3.88 

8 

Pseudorca  crassidens 

Male 

762.9 

3.99 

8 

Delphinapterus  leucas 

Male 

589.7 

3.56 

12 

Delphinapterus  leucas 

Female 

698.5 

3.25 

8 

Delphinapterus  leucas 

Male 

725.7 

4.12 

13 

Globicephala  melaena 

Female 

764.3 

3.99 

15 

Orcinus  orca 

Male 

841.9 

3.81 

1.8 

Orcinus  orca 

Female 

929.9 

4.16 

3 

Orcinus  orca 

Female 

1424.3 

4.46 

16 

Orcinus  orca 

Female 

1995.8 

5.16 

16 

Orcinus  orca 

Female 

2086.5 

5.57 

16 

Orcinus  orca 

Female 

2331.5 

5.25 

10 

Dolphin  Strouhal  numbers  1635 


7.3-11.0  m,  and  water  temperature  ranged  between  12°C  and 
22°C. 

The  animals  normally  swam  about  0.5-1. Om  below  the 
surface  of  the  water  and  1  m  from  the  pool  walls.  Many  of  the 
animals  were  trained  to  perform  ‘fast  swims’  around  the 
circumference  of  the  main  performance  pool,  thereby 
providing  a  larger  range  of  swim  speeds  to  study.  Initial 
training  involved  instructing  the  animals  to  accelerate  quickly 
and  touch  their  rostrum  to  a  boat  pole  held  several  meters 
in  front  of  them.  Through  successive  approximations  of 
positioning  the  boat  pole  further  and  further  away,  the  dolphins 
were  trained  to  swim  fast  upon  command.  When  the  dolphins 
were  not  swimming  at  what  the  trainers  judged  to  be  peak 
performance,  the  boat  pole  was  slapped  on  the  surface  of  the 
water  to  induce  the  dolphins  to  swim  more  rapidly.  Trainers 
frequently  raced  animals  in  pairs  to  provide  additional 
incentive.  A  wide  variety  of  rewards,  including  tactile 
stimulation,  environmental  enrichment  devices  and  food, 
were  given  for  appropriate  behaviors  on  an  intermittent 
reinforcement  schedule.  In  addition  to  performing  in  shows, 
the  dolphins  in  this  study  participated  in  training,  play, 
relationship,  husbandry  and  exercise  sessions  on  a  regular 
basis.  Approximately  18-20  h  of  their  day  consisted  of 
nonstructural  play,  free  and  rest  time. 

Video  analysis 

A  camcorder  (Sony  CCD-TR81  or  Panasonic  DV-510)  was 
used  to  record  swimming  sequences  of  cetaceans  at  a  rate  of 
30  frames  s“*.  The  camcorder  was  positioned  in  front  of  the 
Plexiglas  wall  of  the  pool,  allowing  for  a  clear  view  of  three 
panels.  Swimming  motions  of  the  animals  were  recorded  as 
they  routinely  swam  and  when  they  were  encouraged  to  swim 
at  maximum  speed.  Sequential  body  and  fluke  positions  were 
determined  directly  from  individual  frames  of  videotape  with 
a  Panasonic  AG-7300  video  recorder  and  video  monitor  or 
were  digitized  using  the  Peak  Motus  video  analysis  system 
(version  4.3.1;  Peak  Performance  Technologies,  Englewood, 
CO,  USA).  Kinematic  data  from  video  records  to  calculate 
Strouhal  number  included  mean  swimming  speed  (t/;  m  s“'), 
fluke  oscillation  frequency  (f;  Hz)  and  peak-to-peak  fluke 
amplitude  (A;  m).  A  is  defined  as  the  maximum  vertical 
displacement  of  the  trailing  edge  of  the  flukes.  Only  video 
sequences  in  which  the  animals  appeared  to  be  swimming 
horizontally  and  at  a  constant  speed  were  used. 

At  Sea  Worlds  in  Orlando  and  San  Antonio,  the  animals 
were  marked  with  zinc  oxide  reference  points  on  the  lateral 
aspects  of  the  caudal  peduncle.  Marks  were  separated  by  a 
measured  distance  of  0.1-0.39  m  and  served  as  the  scale  for 
video  analysis.  At  the  San  Diego  site,  the  animals  were  not 
marked.  The  scale  was  determined  from  a  marked  section  of 
the  Plexiglas  panels  of  the  pool.  U  was  determined  by  dividing 
the  length  of  a  marked  section  through  which  the  dolphin  swam 
by  the  time  that  it  took  the  dolphin  to  swim  across  it.  Time  was 
determined  from  the  frame  rate.  Dolphin  swimming  speed 
measurements  obtained  from  the  video  could  be  accurately 
repeated  to  within  a  few  percent.  This  uncertainty  resulted 


from  the  fact  that  time  was  quantified  by  the  frame  rate,  so 
crossing  points  could  be  off  by  a  fraction  of  a  frame  rate.  At 
the  highest  speeds  of  8  m  s“^  the  8  m  run  resulted  in  30  frames 
between  start  and  finish,  of  which  the  last  two  frames  were 
suspect.  This  could  result  in,  at  most,  a  7%  uncertainty  (two 
frames  out  of  30). 

To  assess  if  the  Plexiglas  panels  and  the  recording  position 
affected  the  calculations  of  U,  video  recordings  of  a  cast 
model  of  a  dolphin  dorsal  fin  were  made  as  it  was  moved  along 
the  normal  swimming  trajectory  of  the  animals.  The  difference 
in  distance  between  the  actual  positions  where  the  cast  fin 
crossed  the  reference  marks  and  that  determined  from  the 
video  recordings  was  insignificant.  /  was  calculated  by 
dividing  the  frame  rate  by  the  number  of  frames  comprising 
a  single  complete  oscillation  of  the  tail.  Again,  the  frame  rate 
limited  accuracy  for  determining  the  tail  oscillation  period. 
For  a  relatively  high  frequency  of  3  Hz,  the  period  was 
determined  from  10  frames.  With  a  full  frame  ambiguity  at 
each  end  of  the  oscillation  the  uncertainty  would  be,  at  most, 
20%.  The  peak-to-peak  amplitude  of  the  fluke.  A,  along  with 
a  previously  measured  reference  length  marked  on  the  tank 
wall,  was  measured  directly  on  the  television  monitor  screen. 
The  reference  length  was  recorded  inside  and  outside  the  pool 
to  account  for  refraction  effects.  The  reference  length 
provided  a  means  of  converting  lengths  measured  on  the 
monitor  screen  to  actual  distances  in  meters.  The  spatial 
resolution  of  the  fluke  at  maximum  and  minimum  amplitude 
is,  as  for  the  swim  speed  measurements,  affected  by  the 
framing  rate  of  the  camera.  However,  because  the  fluke  tip  is 
moving  relatively  slowly  at  these  extremes,  the  majority  of  the 
amplitude  uncertainty  resulted  from  resolving  on  the  video 
screen  the  position  of  the  tip  of  the  fluke,  particularly  if  the 
fluke  was  close  to  the  water’s  surface.  When  Strouhal  number 


Table  2.  Regression  equations  for  kinematics 


Species 

Regression  equation 

y2 

Delphinapterus  leucas 

^0.136-l-1.462t//L 

0.861 

n=3 

A/L=0.226-0.034f//L 

0.317 

Globicephala  melaena 

^2.246-0.1 1217/L 

0.318 

n=\ 

A/L=0.086-H0.056t//L 

0.113 

Lagenorhynchus  obliquidens 

^-0.484-h1.25617/L 

0.318 

n=2 

A/L=0.132-i-0.030U/L 

0.077 

Orcinus  area 

7^0.890-i-0.615U/L 

0.529 

n=6 

AIL=QA19+Q.Q2AUIL 

0.066 

Pseudorca  crassidens 

7^0.368-1-1. 056U/L 

0.642 

n=4 

A/L=0.160-i-0.026U/L 

0.095 

Stenella  frontalis 

7^0.408-1-1. 131 U/L 

0.829 

n=\ 

A/L=0.248-0.001f//L 

0.002 

Tursiops  truncatus 

/=0.624-i-0.844U/L 

0.760 

n=\  1 

A/L=0.624-l-0.844t//L 

0.097 

Pooled/ 

/=0.585-i-0.885U/L 

0.802 

Pooled  AIL 

n=28 

A/L=0.184-i-0.015U/L 

0.091 

/  frequency  of  propulsive  cycle;  AIL,  length-specific  peak-to-peak 
amplitude;  LUL,  length-specific  swimming  speed. 
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Fig.  1  (A)  Non-dimensional  fluke-beat  amplitude,  AIL,  as  a  function  of 
length-specific  swimming  speed,  UIL,  where  A  is  the  peak-to-peak  fluke- 
beat  amplitude,  L  is  the  body  length  and  U  is  the  swimming  speed. 
(B)  Non-dimensional  fluke-beat  amplitude,  AIL,  as  a  function  of  fluke- 
beat  frequency  /  (Hz),  where  A  is  the  peak-to-peak  fluke-beat  amplitude 
and  L  is  the  body  length.  Symbols  indicate  particular  cetacean  species: 
Tiirsiops  truncatus  (solid  red  circles),  Pseudorca  crassidens  (open  green 
squares),  Orcinus  area  (solid  blue  triangles),  Globicephala  melaena  (blue 
crosses),  Lagenorhynchus  obliquidens  (solid  blue  diamonds),  Stenella 
frontalis  (open  brown  diamonds)  and  Delphinapterus  leucas  (solid  black 
triangles). 


calculations  were  repeated  for  the  same  recordings  but 
by  independent  observers,  swimming  speed  and  tail 
beat  frequencies  showed  excellent  agreement;  for  tail 
beat  amplitudes,  differences  of  10-20%  were  not 
uncommon.  This  uncertainty  resulted  from 
insufficient  screen  resolution,  framing  rate  and  the 
proximity  of  the  tail  to  the  water  surface.  Overall  St 
uncertainty  was  estimated  to  be  -20%.  This  is  6%  less 
than  the  St  uncertainty  calculated  for  the  worse  case 
scenario  by  propagating  the  independent  uncertainties 
estimated  for  high  speeds  (7%),  high  frequencies 
(20%)  and  a  fluke  amplitude  uncertainty  of  15%. 

To  adjust  for  size  differences  between  species,  data 
were  analyzed  with  respect  to  length-specific  velocity 
{UIL)  and  length-specific  amplitude  {AIL).  In  some 
analyses,  /  was  non-dimensionalized  by  dividing 
frequency  by  UIL.  Means  were  calculated  for  values 
that  did  not  vary  with  L  or  UIL.  Variation  about  means 
was  expressed  as  +1  S.D.  Linear  relationships  were 
estimated  by  least-squares  regression  (Microsoft 
Excel).  Differences  in  slopes  of  the  regressions  were 
analyzed  by  analysis  of  covariance  (Zar,  1984).  Means 
of  species  were  compared  using  analysis  of  variance 
(ANOVA;  Statistica  Version  4.1,  StatSoft).  Results 
were  considered  significant  at  the  a=0.05  level. 

Results 

Kinematic  data 

A  total  of  267  swimming  sequences  were  used  for 
kinematic  analysis  {S.  frontalis,  n=13;  L.  obliquidens, 
n=17;  T.  truncatus,  n=107;  P.  crassidens,  n=69;  G. 
melaena,  n=l2;  O.  orca,  n=30;  D.  leucas,  n=\9).  The 
animals  maintained  continuous  propulsive  motions  by 
vertical  oscillations  of  the  flukes,  as  has  been  described 
previously  (Fish,  1993,  1998;  Rohr  et  ah,  2002). 
The  fastest  mean  swimming  speeds  were 
t/=6.42+0.41  m  s-'  and  t//L=2.90+0.19  L  s”'  for  L. 
obliquidens  and  the  slowest  mean  speeds  were 
(7=2.38+0.74  m  s“'  and  (//L=0. 68+0.22  L  s“*  for  D. 
leucas. 


Fig.  2.  Fluke-beat  frequency,  /  (Hz),  as  a 
function  of  length-specific  swimming  speed, 
UIL,  where  L  is  the  body  length  and  U  is  the 
swimming  speed.  Regression  equation  for  all 
the  cetacean  data  is:  y=0.89x-l-0.59  (F=0.8, 
N=261).  Symbols  indicate  particular  cetacean 
species:  Tursiops  truncatus  (solid  red 
circles),  Pseudorca  crassidens  (open  green 
squares),  Orcinus  orca  (solid  blue  triangles), 
Globicephala  melaena  (blue  crosses), 
Lagenorhynchus  obliquidens  (solid  blue 
diamonds),  Stenella  frontalis  (open  brown 
diamonds)  and  Delphinapterus  leucas  (solid 
black  triangles). 
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With  no  statistically  significant  correlation,  AIL  was  found 
to  be  relatively  insensitive  to  both  UIL  (Fig.  1  A)  and/(Fig.  IB) 
for  all  species.  The  mean  value  of  AIL  for  all  odontocetes  was 
0.21+0.03  (n=267).  Mean  AIL  ranged  from  0.25+0.02  for  S. 
frontalis  to  0.17+0.02  for  G.  melaena  with  89%  of  the  data 
residing  between  0.15  and  0.25.  ANOVA  showed  that  there 
was  a  signihcant  difference  for  AIL  among  species  (7’<0.001; 
7^=9. 76;  d.f.=6,  260).  Aggregating  all  the  odontocete  data 
(n=267),  /  was  found  to  increase  linearly  with  increasing  UIL 
as  /=0.89(t//L)+0.59  (r^=0.S;  Fig.  2).  A  positive  linear 
relationship  between/and  UIL  is  similar  to  results  reported  for 
cetaceans,  hsh  and  other  marine  mammals  (Bainbridge,  1958; 
Hunter  and  Zweifel,  1971;  Webb  and  Kostecki,  1984; 
Feldkamp,  1987;  Fish  et  ah,  1988;  Scharold  et  ah, 1989). 
Regression  equation  for/and  AIL  with  respect  to  UIL  for  each 
species  is  provided  in  Table  2.  The  negative  slope  in  the 
regression  equation  for /by  G.  melaena  is  due  to  the  limited 
speed  range. 

Strouhal  data 

The  computed  Strouhal  number  showed  little  dependence  on 
body  length  or  swim  speed  for  the  delphinid  species  (Fig.  3). 
Aggregating  animals  for  each  species  (Fig.  4),  mean  St  values 
generally  reside  near  the  lower  boundary  of  the  0.25-0.35 
range  (Table  3)  predicted  by  Triantafyllou  et  al.  (1991,  1993) 


Table  3.  Mean  Strouhal  values 


Sample  size 

Species 

Mean  +  s.D. 

n=107 

Tursiops  truncatus 

0.26±0.05 

n=69 

Pseudorca  crassidens 

0.26±0.05 

n=30 

Orcinus  orca 

0.28+0.05 

11=19 

Delphinapterus  leucas 

0.35+0.1 

11=17 

Lagenorhynchus  obliquidens 

0.24+0.03 

11=13 

Stenella  frontalis 

0.33+0.03 

11=12 

Globicephala  melaena 

0.24+0.02 
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Fig.  3  (A)  Strouhal  number,  St,  of  odontocete  cetaceans  as  a  function  of  body  length,  L. 
(B)  Strouhal  number,  St,  of  odontocete  cetaceans  as  a  function  of  swimming  speed,  U. 
Symbols  indicate  particular  species:  Tursiops  truncatus  (solid  red  circles),  Pseudorca 
crassidens  (open  green  square),  Orcinus  area  (solid  blue  triangles),  Globicephala 
melaena  (blue  crosses),  Lagenorhynchus  obliquidens  (solid  blue  diamonds)  and 
Stenella  frontalis  (open  brown  diamonds). 


for  peak  propulsive  efficiency.  Excluding  D.  leucas,  the  mean 
St  for  the  delphinids  was  0.26+0.05  (n=248).  The  predicted 
0.25-0.35  St  range  captured  55%  of  the  delphinid  St  data 
(Fig.  5),  whereas  the  range  from  0.2  to  0.3  contained  74%.  For 
an  incremental  St  range  of  0.05,  the  majority  of  the  data  were 
found  between  0.225  and  0.275  (44%).  D.  leucas  had  a  mean 
St  of  0.35+0.10  (n=19),  which  was  conspicuously  higher  than 
most  of  the  St  values  for  the  delphinids. 

Propulsive  efficiencies,  which  were  previously  reported  by 
Fish  (1998),  are  plotted  as  a  function  of  St  in  Fig.  6.  For  P. 
crassidens,  O.  orca  and  L.  truncatus,  propulsive  efficiencies 
were  found  to  broadly  peak  at  about  0.90,  0.87  and  0.85, 
respectively,  over  a  relatively  narrow  range  of  St  (0.23-0.28). 
Outside  this  St  range,  where  measurements  exist,  efficiencies 
drop  off  rapidly.  The  St  range  favored  by  P. 
crassidens,  O.  orca  and  T.  truncatus  was 
within  this  same  range,  0.225<5'f<0.30 
(Fig.  7A-C).  The  efficiency  data  for  D.  leucas 
were  lower  (0.83)  and  exhibit  a  conspicuously 
broader  peak  at  5'f=0.25-0.40.  The  distribution 
4  of  D.  leucas  St  was  relatively  flat,  with  a 

narrow  peak  occurring  at  5'f=0.425-0.45 
(Fig.  7D). 


Discussion 

Cetaceans  swim  by  oscillatory  heaving  and 
pitching  of  the  caudal  flukes,  which  act  as  a 
hydrofoil  (Lighthill,  1969;  Webb,  1975;  Fish 
and  Hui,  1991;  Fish,  1993,  1998).  The 
oscillating  movements  of  a  hydrofoil  result  in 
unsteady  shedding  of  vorticity  from  the  trailing 
edge  (Anderson  et  al.,  1998).  The  pattern  and 
spin  of  the  staggered  array  of  vortices  generate 
a  jet  flow,  which  produces  thrust  to  overcome 
the  drag  on  the  body.  Triantafyllou  et  al.  (1991, 
1993,  2002)  considered  the  jet  to  be 
convectively  unstable,  acting  as  a  tuned 
ampliher  with  a  narrow  range  of  frequencies  of 
maximum  ampliheation  (i.e.  maximum  thrust 
production).  The  pattern  and  periodicity  of 
vortices  shed  into  the  wake,  therefore, 
determine  the  optimal  thrust  production  for 
maximum  efficiency.  The  arrangement  of 
vortices  for  maximum  efficiency  is  a  reverse 
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Fig.  4.  Strouhal  number  (St)  range  (+  s.D.)  for  cetaceans. 
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Fig.  5.  Flistogram  of  percentage  of  delphinid  cetacean  Strouhal  number  (St)  data 
(n=248);  incremental  St  range  is  0. 1,  incremental  step  is  0.25. 


Fig.  6.  Propulsion  efficiency  of  cetaceans  as  a  function  of  Strouhal  number  (St). 
Colors  indicate  particular  species:  Tursiops  truncatus  (red),  Pseudorca  crassidens 
(green),  Orcinus  area  (blue)  and  Delphinapterus  leucas  (black). 


Karmen  street  (Triantafyllou  et  al.,  2002).  The 
wake  dynamics  are  dominated  by  the  non- 
dimensional  Strouhal  number,  St,  in  which  the 
distance  between  vortices  and  their  rate  of 
formation  co-vary  with  speed  (Vogel,  1994; 
Triantafyllou  and  Triantafyllou,  1995; 
Triantafyllou  et  al.,  2002).  Experimental  studies 
of  heaving  and  pitching  foils  have  found  that  the 
structure  of  the  vortex  wake  changes  with  St  (von 
Ellenrieder  et  al.,  2003;  Taylor  et  al.,  2003), 
maximum  thrust  occurred  between  0.25  and  0.4 
(Triantafyllou  et  al.,  1993;  Anderson  et  al.,  1998) 
and  maximum  efficiency  was  within  the  range  of 
0.25-0.4  (Triantafyllou  et  al.,  1991,  1993; 
Anderson  et  al.,  1998).  Efficiencies  for  flapping 
foil  experiments  have  been  reported  to  peak 
below  the  optimum  0.25-0.35  St  range  predicted. 

The  data  for  all  delphinids  (u=248)  showed 
little  dependence  on  St  over  a  range  of  swim 
speeds  from  about  2-8  m  s“*  (Fig.  3B).  It  has 
been  hypothesized  for  cruising  flight  and 
swimming  that  St  would  be  ‘tuned’  for  high 
propulsive  efficiency  (Triantafyllou  et  al.,  1991, 
1993;  Taylor  et  al.,  2003).  Cruising  speeds  for 
the  cetaceans  have  been  reported  from 
~l-5  m  s“'  (see  Fish,  1998).  For  the  present  St 
data  (Fig.  5),  a  conspicuous  peak  was  not 
apparent  at  cruising  speeds  (Fig.  3B)  or  where 
maximum  propulsive  efficiency  was  predicted 
(St=0.25-0.35).  Moreover,  the  St  data  were  not 
most  concentrated  in  the  predicted  range. 
Whereas  55%  of  the  data  fell  within  the 
predicted  range  of  5'f=0.25-0.35,  74%  of  all  the 
St  data  occurred  between  0.2  and  0.3.  Some 
of  the  scatter  in  the  St  data  is  a  result  of 
measurement  uncertainty.  Wolfgang  et  al. 
(1999)  have  reported  St  uncertainties  of  -30% 
for  studies  with  fish  (Danio  malabaricus). 
However,  a  large  part  is  presumably  due  to 
natural  variation  of  the  animal’s  swimming 
motion  (Rosen,  1959;  Wolfgang  et  al.,  1999). 
Kayan  and  Pyatetskiy  (1977)  reported  a 
dependence  of  St  on  acceleration  for  captive  T. 
truncatus,  with  St  increasing  with  increasing 
acceleration.  Taylor  et  al.  (2003)  similarly  found 
that,  for  birds,  St  was  significantly  higher  for 
intermittent  as  opposed  to  direct  flight.  Although 
data  from  the  present  video  analysis  were  limited 
to  steady  swimming  speeds,  effects  due  to  small 
accelerations  were  possible.  A  dependence  of 
St  on  acceleration  may  partly  explain  the 
difference  in  St  values  for  L.  obliquidens 
between  the  present  data  (5'f=0.24+0.03,  n=\l) 
and  those  inferred  from  traces  of  an  accelerating 
animal  (5'f=0.30,  n=l;  Triantafyllou  et  al.,  1991, 
1993). 
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Fig.  7.  Histogram  of  percentage  of  (A)  Pseudorca  crassidens  (n=69),  (B) 
Orcinus  area  (n=30),  (C)  Tursiops  truncatus  {n=\Ql)  and  (D)  Delphinapterus 
leucas  (n=19)  Strouhal  number  (St)  data. 


One  cannot  be  certain  that  the  mode  of  steady 
swimming  in  captivity  for  relatively  short  durations 
and  near  the  water  surface  is  similar  to  that 
employed  for  long  durations  in  the  wild.  However, 
the  high-speed  swimming  capability  of  regularly 
exercised  captive  and  free-ranging  dolphins  is 
generally  similar  (Rohr  et  al.,  2002).  Pershin  (1988) 
reports  an  St  value  of  0.37  for  a  free-ranging 
dolphin,  which  he  refers  to  as  a  common  dolphin. 

Pershin  (1988)  makes  no  reference  to  whether  the 
animal  was  accelerating  or  how  the  recordings  were 
made.  Unlike  the  captive  dolphins  in  the  present 
study,  this  free-ranging  dolphin  was  not  swimming 
near  the  surface.  Except  for  D.  leucas  and  S. 
frontalis,  an  St  value  of  0.37  is  conspicuously  higher 
than  the  mean  values  reported  here.  It  is  not  known 
if  this  disparity  reflects  differences  between  species, 
captive  and  free-ranging  animals,  steady  or 
accelerated  swimming  or  different  depths  beneath 
the  surface. 

A  hydromechanical  model  of  lunate-tail 
propulsion  based  on  three-dimensional  unsteady 
wing  theory  with  continuous  loading  (Chopra  and 
Kambe,  1977;  Yates,  1983)  was  used  by  Fish  (1998) 
to  calculate  efficiency  for  O.  area,  T.  truncatus,  P. 
crassidens  and  D.  leucas.  Here,  efficiency  refers  to 
the  Froude  efficiency  defined  as  the  mean  rate  of 
mechanical  work  derived  from  mean  thrust,  divided 
by  all  the  work  that  the  animal  is  performing  while 
swimming  (Chopra  and  Kambe,  1977).  The 
calculated  efficiencies  are  similar  to  values  reported 
previously  using  different  hydrodynamic  models 
(Fish,  1998).  Maximum  efficiencies  of  85-90% 
generally  occurred  at  typical  cruising  speeds  (Fish, 

1998).  Similarly,  the  minimum  cost  of  transport 
coincides  with  the  range  of  cruising  speeds  in  T. 
truncatus  (Williams  et  al.,  1992;  Yazdi  et  al.,  1999). 

Efficiency  values  as  high  as  85%  have  been 
measured  for  advanced  torpedo  propellers  that  take 
advantage  of  the  boundary  layer  wake  (Fang  and 
Daybell,  1963). 

Dolphin  efficiencies  were  significantly  higher 
than  values  measured  for  conventional  small  rigid 
propellers,  which  are  typically  no  more  that  40% 
(Triantafyllou  and  Triantafyllou,  1995).  The  higher 
efficiencies  of  dolphins  are  probably  due  to  some 
degree  of  St  ‘tuning’  (in  0.2  St  increments,  94%  of 
the  St  data  were  between  0.2  and  0.4),  oscillatory 
fluke  motions  and  greater  flexibility  of  the  flukes. 
Oscillating  foils  can  produce  high  lift  coefficients 
and  efficiencies  by  vorticity  control  (Gopalkrishnan  et  al., 
1994;  Anderson  et  al.,  1998;  Ramamurti  et  al.,  2002;  Read  et 
al.,  2003).  The  unsteady  effects  can  increase  the  lift  and  permit 
the  foil  to  function  at  high  angles  of  attack  without  stalling. 
Cetacean  flukes  are  lateral  extensions  of  the  tail  and  are  a 
composite  of  flexible  elements,  including  dense  fibrous  tissue. 


ligaments  and  blubber  (Felts,  1966).  The  only  rigid 
components  in  the  cetacean  tail  are  a  series  of  short  caudal 
vertebrae,  which  support  the  flukes  axially  (Rommel,  1990). 
The  white-sided  dolphin  (Lagenorhynchus  acutus)  shows  35% 
and  13%  chordwise  and  spanwise  deflections,  respectively 
(Curren  et  al.,  1994).  Chordwise  flexibility  of  an  oscillating  foil 
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has  been  demonstrated  experimentally  and  theoretically  to 
increase  efficiency  by  up  to  36%  with  only  a  small  reduction 
in  thrust  compared  with  a  rigid  foil  (Katz  and  Weihs,  1978, 
1979;  Bose,  1995;  Prempraneerach  et  al.,  2003).  Heaving  and 
pitching  motions  of  flexible  foil  produced  a  maximum 
efficiency  of  0.87  at  St=Q3,  and  an  efficiency  of  0.8  was 
achieved  at  various  combinations  of  angles  of  attack  around 
15°  and  St  ranging  from  0.17  to  0.35  (Prempraneerach  et  al., 
2003). 

Prior  to  the  present  investigation  there  have  been  no  studies, 
for  either  cetaceans  or  fish,  that  addressed  whether  swimming 
efficiency  occurs  within  the  predicted  0.25-0.35  range  of 
Triantafyllou  et  al.  (1991,  1993)  (Bandyopadhyay  et  al.,  2000). 
Cetacean  St  values  occurred  in  the  present  study  most 
frequently  in  the  range  0.225<5f<0.275  (Figs  5,  7A-C).  Peak 
Froude  efficiencies  were  found  in  this  same  range  (Fig.  6), 
which  straddles  the  lower  boundary  of  the  0.25-0.35  St  range 
predicted.  Although  there  are  no  cetacean  or  fish  efficiency 
measurements  to  compare  with,  there  are  relevant  data  for  two- 
dimensional  foils.  Bandyopadhyay  (2002)  reported  the  peak 
efficiency  for  a  pair  of  tail  flapping  foils  to  be  below  the 
0.25-0.35  St  range.  The  efficiency  for  oscillating  foils  reported 
by  Triantafyllou  et  al.  (1993),  although  maximal  in  the  St  range 
between  0.2  and  0.35,  is  practically  flat.  The  efficiency  versus 
St  data  of  Anderson  et  al.  (1998;  Fig.  5)  for  a  two-dimensional 
foil  clearly  do  not  show  peak  efficiencies  tuned  to  an  St  range 
between  0.25  and  0.35.  The  lack  of  closer  agreement  with  the 
present  data  is  not  surprising  given  that  the  foil  experiments 
could  not  (and  were  not  assumed  to)  capture  the  full  flow  field 
dynamics  of  swimming  animals. 

St  is  often  expressed  as  the  product  of  non-dimensional 
peak-to-peak  fluke  amplitude  (A/L)  and  non-dimensional 


Fig.  8.  Non-dimensional  fluke-beat  amplitude,  A/L,  as  a  function  of  non- 
dimensional  fluke-beat  frequency,  f/(U/L),  where  A  is  the  peak-to-peak 
fluke  amplitude,  L  is  the  body  length,  /  is  the  fluke-beat  frequency  and  U  is 
the  swimming  speed.  Contours  of  constant  Strouhal  number,  St,  are 
included.  Symbols  indicate  particular  cetacean  species:  Tursiops  truncatus 
(solid  red  circles),  Pseudorca  crassidens  (open  green  squares),  Orcinus 
orca  (solid  blue  triangles),  Globicephala  melaena  (blue  crosses), 
Lagenorhynchiis  obliquidens  (solid  blue  diamonds),  Stenella  frontalis  (open 
brown  diamonds)  and  Delphinapterus  leucas  (solid  black  triangles). 


frequency  [/7(t//L)].  Plotting  the  data  in  these  coordinates  and 
including  contours  of  constant  St  values  provides  an  additional 
perspective  on  the  parameters  controlling  the  range  of  St 
encountered  (Fig.  8).  Regardless  of  St  value,  89%  of  all  the  A/L 
data  («=267)  falls  within  a  range  of  0.15-0.25.  Cetaceans  are 
known  to  maintain  a  nearly  constant  A/L  of  0.20  with  respect 
to  U  (Kayan  and  Pyatetskiy,  1977;  Fish,  1998;  Fish  et  al., 
2003).  A  value  of  0.20  for  AIL  is  also  typical  for  other  animals 
that  swim  by  movements  of  the  body  and  tail  (Bainbridge, 
1958;  Webb,  1975;  Videler,  1993;  Fish,  1998;  Schultz  and 
Webb,  2002;  Fish  et  al.,  2003). 

Swim  speed  is  controlled  through  frequency  modulation,  as 
is  common  for  fish  and  other  marine  mammals  that  swim  by 
oscillations  of  high  aspect  ratio  hydrofoils  (Bainbridge,  1958; 
Hunter  and  Zweifel,  1971;  Webb  and  Kostecki,  1984;  Feldkamp, 
1987;  Fish  et  al.,  1988;  Scharold  et  al.,  1989).  Regardless  of  St 
value,  90%  of  fl{UIL)  data  falls  within  a  range  of  1-2.  Note,  if 
A/L=0. 15-0.25  and//(t//L)=l-2,  the  corresponding  St  range  is 
0.15-0.50.  Presumably  the  ‘boundaries’  imposed  by  AIL  and 
fl{UIL)  on  St  are  manifestations  of  additional  morphological  and 
hydrodynamic  constraints  imposed  on  the  animal.  Optimal  St 
values  for  pitching,  heaving  and  flapping  foils  have  also  been 
found  to  depend  on  other  kinematic  parameters,  including  angle 
of  attack,  amplitude-to-chord  ratio  and  phase  of  motion 
(Anderson  et  al.,  1998;  Wang,  2000;  Bandyopadhyay  et  al., 
2000;  Read  et  al.,  2003). 

Representative  St  has  often  been  estimated  by  simply 
multiplying  the  slope  of  frequency  versus  speed  over  length 
by  the  mean  non-dimensional  peak-to-peak  fluke  amplitude 
(Triantafyllou  et  al.,  1993;  Bandyopadhyay,  2002).  This 
assumes  the  /  versus  LUL  data  passes  through  the  origin.  The 
slope  of  /  versus  LUL  throughout  the  present  data  is  0.89 
(Fig.  2).  The  mean  AIL  is  0.21.  This  product  results  in 
an  estimate  of  0.19  for  St  for  odontocete  cetaceans. 
However,  if  the  product  of  fl(UIL)  and  AIL  is  first 
calculated  for  each  observation  and  then  averaged, 
St=0.21 .  A  30%  disparity  results  because  the  y-intercept 
(0.58;  see  Fig.  2)  of  the  least  squares  estimate  of /versus 
LUL  was  initially  ignored. 

Whereas  91%  of  all  the  cetacean  data  fall  between 
St=0.2  and  0.4,  a  significant  fraction  (37%)  of  the  D. 
leucas  data  conspicuously  fell  outside  this  range.  Mean 
St  for  D.  leucas  is  0.35+0.10  compared  with  0.27+0.06 
for  the  other  cetaceans  examined  in  this  study.  It  has 
been  noted  that  the  swimming  performance  and 
efficiency  of  D.  leucas  differ  considerably  from  those  of 
other  cetaceans,  which  is  consistent  with  its  general 
body  contour  and  low  aspect  ratio  flukes  (Fish,  1998). 
D.  leucas  generally  exhibits  the  poorest  swimming 
performance  of  cetaceans.  This  species  feeds  on  slow- 
moving  prey,  including  crustaceans  and  annelids 
(Brodie,  1989).  The  mean  St  for  S.  frontalis  was  also 
high  at  0.33+0.03.  Unlike  D.  leucas,  S.  frontalis  is  a  fast 
swimmer  (Fish  and  Rohr,  1999).  The  individual  S. 
frontalis  examined  in  this  study  had  a  mean  AIL  of 
0.25+0.02,  which  was  the  highest  of  any  species  tested. 


The  data  presented  in  this  study  significantly  expand  the 
previously  used  Strouhal  data  for  cetaceans  (Triantafyllou  et 
ah,  1991,  1993;  Triantafyllou  and  Triantafyllou,  1995;  Taylor 
et  al.,  2003)  by  greatly  increasing  the  number  of  observations 
(more  than  100-fold),  species  (7-fold)  and  range  of  Reynolds 
number  (10-fold)  and,  for  the  first  time,  provide  a  direct 
comparison  between  measured  Strouhal  number  and 
swimming  efficiency.  The  present  data  show  that  over  the 
range  of  swim  speeds  observed,  2-8  m  s“*,  cetaceans  swim  at 
St  values  between  0.2  and  0.4,  preferring  a  range  of  0.2-0. 3 
where  maximum  efficiencies  occur.  The  strong  relationship 
between  maximum  propulsive  efficiency  and  St  continues  to 
support  the  premise  (Triantafyllou  et  al.,  1991,  1993)  that 
vorticity  control  associated  with  fluke  kinematics  is  an 
important  attribute  of  cetacean  swimming  performance. 
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